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Hybrid Monte Carlo simulations of the pure compact U(l) gauge theory are performed with the Tsallis weight. 
The simulations show that the use of the Tsallis weight enhances the tunneling rate between metastable states. 



1. INTRODUCTION 

The four dimensional pure compact U(l) gauge 
theory is known to posses a phase transition (PT) 
separating a confined phase and a Coulomb one 
and the determination of the order of PT is of 
great importance. In the pure compact U(l) 
gauge theory, however, the determination of the 
order of PT is very difficult since on large lattices 
the standard updating algorithms like Metropo- 
lis, heat-bath, hybrid Monte Carlo (HMC)||l|, fail 
to generate enough tunneling between metastable 
states. Many efforts have been done to clarify the 
order of the PT|^. In spite of such efforts the or- 
der of the PT is still controversial. 

A promising algorithm to overcome this dif- 
ficulty is the multicanonical algorithm[^ which 
uses a multicanonical weight in stead of the Boltz- 
mann one and can enhance the tunneling rate 
between metastable states. Although the multi- 
canonical algorithm works effectively, one disad- 
vantage of the multicanonical algorithm might be 
that a multicanonical weight used in the simula- 
tions is not known a priori and it must be deter- 
mined before the simulations. Usually the weight 
is estimated from a short run of the standard 
updating algorithm. However this estimation is 
more difficult as a lattice size becomes larger. 

Recently in several fields where the usual 
Monte Carlo technique does not work effec- 
tively, simulations were performed with the Tsal- 
lis weight first introduced by TsaUis||]. The 
Tsallis weight controlled by a parameter q can be 
easily defined, and the usual Boltzmann weight is 
given by taking the limit q I. The results have 
showed that the use of the Tsallis weight might be 
more advantageous than that of the Boltzmann 



one. 

Inspired by such studies, here we apply the 
Tsallis weight for the hybrid Monte Carlo sim- 
ulation of the pure compact U(l) lattice gauge 
theory and investigate whether the Tsallis weight 
solves the difficulty in determining the order of 
the PT. 

2. GENERALIZED ENSEMBLE AND 
TSALLIS WEIGHT 

The partition function of a system with the ac- 
tion S is given by 
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Introducing a function G, the partition function 
can be rewritten as 



Z ^ I eiip{~pS[U])/G-GdU 
exp(-/3S') 
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where Zq is a generalized partition function: 

Zg = j GdU, (3) 

and < O >G stands for the expectation value of 
an observable O in the ensemble generated with 
Zg- The expectation value of O in the canonical 
ensemble generated with cq.(Q) can be obtained 
in terms of the generalized ensemble as 
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Previously it was suggested that in some cases 
Tsallis weight works better than the Boltzmann 
onelH.H. In the TsaUis formulation, 



G = [1 - {1 - q)f3S] — , 

or in the familiar exponential form, 

G = exp(-/3S'T) = Wt, 



with 



HI - (1 - q)(3S], 



(5) 



(6) 
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where S = S—So- 5*0 is a certain constant shifting 
the origin of the action. Here we call Wt " Tsallis 
weight" . The Boltzmann weight can be obtained 
by taking the limit g — > 1 in the Tsallis weight 
Wt- 

3. Hybrid Monte Carlo with Tsallis weight 

The action of pure compact U(l) lattice gauge 
theory is given by 



E 

X.fJ,>V 



5= > ' [l-cos(0^,(x))], 



(8) 



where 9^^, (x) is the sum of link angles contribut- 
ing a plaquette at a site x on the fiv plane. The 
Hamiltonian with the Tsallis weight, used for 
HMC simulations, is defined by 



(9) 



Here note that St is not always well-defined. 
Namely 1 — (1 — q)(3S can be negative. However, 
provided that q is close to one, 1 — (1 — q)pS is 
made positive. 

4. RESULTS 

As an exploratory study we tested the above al- 
gorithm varying q on an 8^ lattice in the vicinity 
of the PT ( at /3 = 1.00737 ). We performed HMC 
simulations^] with eq.(^). 5*0 ( = 0.37 ) was taken 
to be the average value of the action S. In order to 
make comparisons, we also performed HMC sim- 
ulations with the Boltzmann weight. Typically 
we simulated 10^ trajectories for each q. 
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Figure 1. (top): Typical history of action S for 
q = 1 { Boltzmann weight ). (bottom): St for 
q=0.99997. 

We observed more tunneling for q < 1 than 
for q = 1 (the Boltzmann weight). Fig.l shows 
representative history of actions of q=l and 
0.999997 and we see more tunneling for the case 
of q=0.999997. 

Fig. 2 shows action densities of q=l, 0.999997, 
0.999995 and 0.99999. The height of the valley 
between two peaks increases as q decreases which 
indicates that the tunneling rate increases as q 
decreases. For too small q, however, values of 
the action will locate near Sq and in this case 
there may be a difficulty in reweighting. Fig. 3 
and 4 show that the action densities reweighted 
to q=l and compare those with the Boltzmann 
action density ( q = 1 ). The reweight densi- 
ties reproduce well the Boltzmann action density. 
However, for q = 0.99999 ( the smallest q here ) 
the reweighted action density has large fluctua- 
tions at S far from Sq, which may indicate that 
one needs caution when using the Tsallis weight. 
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It is also possible to take q > 1 in the simula- 
tions. For g > 1, however, the desired behavior, 
i.e. more tunneling, was never observed for our 
case. 
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Figure 2. The action densities N{S) for q=l, 
0.99997,0.999995 and 0.99999 
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Figure 3. Comparison between the action den- 
sity of q = 1 and the reweighted action density of 
q = 0.999997. 

5. DISCUSSION 

The exploratory study on an 8"* lattice showed 
that the Tsallis weight enhances the tunneling 
rate. The lattice size used here is small. For pre- 
cise determination of the order of the PT, larger 
lattices are needed. Therefore it is important to 
confirm that this preferred feature of the Tsallis 
weight is preserved for larger lattices, possibly up 
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Figure 4. Comparison between the action den- 
sity oi q = I and the reweighted action density of 
q = 0.99999. 

to 18^ lattice size or so. Simulations on larger 
lattices are in progress. 
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